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To explore the electronic structure of the first aromatic superconductor, potassium-doped 
solid picene which has been recently discovered by Mitsuhashi et al. with the transition tem- 
peratures T c — 7 — 20 K, we have obtained a first-principles electronic structure of solid picene 
as a first step toward the elucidation of the mechanism of the superconductivity. The un- 
doped crystal is found to have four conduction bands, which are characterized in terms of the 
maximally localized Wannier orbitals. We have revealed how the band structure reflects the 
stacked arrangement of molecular orbitals for both undoped and doped (Kapicene) cases, where 
the bands are not rigid. The Fermi surface for Kapicene is a curious composite of a warped 
two-dimensional surface and a three-dimensional one. 
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Introduction — While there are mounting interests to- 
ward superconductivity in more versatile classes in recent 
years, fascination continues with superconductivity in n- 
electron (organic) systems. This has a long history, dat- 
ing back to superconductivity in graphite intercalation 
compounds (GIC's), first discovered 1 in 1965 in alkal- 
metal doped ACs (A = K, Rb, Cs) with transition tem- 
peratures T c < 1 K, which was recently raised to 11.5 K 
in CaCg. 2 The superconductivity in organic metals was 
first found in TMTSF compounds 3 in 1980 with T c = 0.9 
K under 12 kbar. Then comes the BEDT-TTF family, 4 
for which T c = 14.2 K under 82 kbar was found in 2003, 
the highest T c to date among the organic superconduc- 
tors. Another family is the doped fullerides, discovered 5 
in 1991 for K3C60 with T c = 18 K at ambient pressure, 
followed by Cs2RbCgo with T c — 33 K at ambient pres- 
sure 6 and CS3C60 with T c — 40 K under 15 kbar, 7 which 
is the highest T c among carbon-based materials. How- 
ever, a new class of superconducting 7r-electron materials 
has not been found over a decade. 

So a recent discovery of superconductivity by Mit- 
suhashi et al. 8 in a solid of aromatic molecule picene, 
with T c = 7 and 20 K when doped with potassium, 
is a long-awaited breakthrough. The discovery is semi- 
nal in that this is the first aromatic superconductor, so 
that a comparison with other classes of 7r-electron su- 
perconductors should be interesting. A picene molecule 
consists of five benzene rings joined in an armchair man- 
ner with a higher chemical stability than pentacene with 
a straight arrangement. The solid picene has in fact a 
much larger band gap (~ 3.3 eV) 9 than that of solid 
pentacene (~ 1.8 eV), which is why the material origi- 
nally attracted attention as an organic semiconductor in 
e.g. organic FET's. 9, 10 Picene has been known to crystal- 
lize in orthorhombic and monoclinic forms, 11 and Kubo- 
zono's group discovered superconductivity in potassium- 
doped monoclinic crystalline picene. 8 The crystal com- 



prises layers stacked in c direction, where each layer (in 
ab plane) has picene molecules arranged in a herring- 
bone structure. The measured lattice parameters of the 
crystalline picene are a — 8.480,6 = 6.154, c = 13.515 
in A with (3 — 90. 46°, 11 where a primitive cell contains 
two incquivalcnt molecules A and B whose positions are 
related with a space group P2\. In 8 Kapicene with the 
concentration x of K atoms per picene is found to become 
a superconductor for x ~ 3, and the Meissner effect has 
been detected with H C 2 > 10 4 Oe suggested. 

Theoretically, the first thing we have to do is first- 
principles electronic structure calculations, which we 
have done in this Letter, with the density functional the- 
ory (DFT) for both pristine and doped solids of picene as 
a first step to elucidate the mechanism of superconduc- 
tivity. We have identified the main characters of the con- 
duction bands in terms of the relevant molecular orbitals, 
and we have also constructed a downfolded four-band 
tight-binding model with maximally localized Wannier 
orbtials. From these we shall show how the band struc- 
ture reflects the stacked arrangement of molecular or- 
bitals for both undoped picene and Kapicene. The Fermi 
surface for Epicene has turned out to be a curious com- 
posite of a warped two-dimensional surface and a three- 
dimensional one, where both the molecular arrangements 
modified by doping and hybridization with alkali-metal 
atoms are relevant. 

First-principles bands — We have first optimized the 
atomic positions including hydrogen atoms with the lat- 
tice parameters fixed at the measured values. 11 For the 
electronic structure calculation for the optimized struc- 
ture we adopt DFT with the projector augmented-wave 
(PAW) method 12 with the Quantum MAterials Simu- 
lator (QMAS) package 13 in the local-density approxi- 
mation (LDA). The pseudo Bloch wave functions are 
expanded by plane waves with an energy cutoff of 40 
Ry and 4 3 fc-points. The angle between the planes of 
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the inequivalent molecules in the optimized geometry is 
61°, close to the measured 58°. Each molecule is not ex- 
actly planar due to the intramolecular overcrowding, as 
seen in the experiment, 11 with the calculated center-of- 
mass vectors between two inequivalent, nearest- neighbor 
molecules re— ta = 0.485a±0.5b. In the calculated elec- 
tronic band structure in Fig.l (a), we first notice that 
the dispersion along the c* axis is small compared to 
those along the a*b* plane, reflecting the layered struc- 
ture. The lowest 102 bands are fully occupied, with the 
valence band top occurring at c*/2. The conduction band 
bottom is also located there, giving a band gap of 2.36 
eV, while the experimental value 9 is 3.3 eV. Such an un- 
derestimation is generally seen in LDA calculations. As 
seen in Fig.l (a), the band dispersions along kf,* = ±1/2 
are doubly degenerate, unlike the triclinic crystalline 
tetracene and pentacene but like the monoclinic crys- 
talline naphthalene and anthracene. 14 The electron den- 
sity, displayed in Fig.l (c), spreads over each molecules 
basically uniformly, with no significant dimerizaton seen. 

For comparison we have also calculated the electronic 
structures of an isolated picene molecule with the molec- 
ular geometries fixed at those optimized in the crystal to 
obtain the HOMO and LUMO wave functions, where the 
LUMO-HOMO energy gap is 2.96 eV (against 1.12 eV 
calculated for pentacene molecule). If we regard picene 
and pentacene as fragments of graphene with armchair 
and zigzag edges, respectively, the higher stability (i.e., 
larger gap) in picene is consistent with the wisdom ob- 
tained in graphene, where armchairs are much stabler. 

If we go back to the crystal, valence bands are the 
99th (VB4) - 102nd (VB1) bands in an energy region of 
width 0.77 eV separated from the lower ones by 0.26 eV. 
Looking at the Bloch wave functions for VB1 and VB2 
at r, we can identify that they have mainly the charac- 
ter of the HOMO's of picene molecule. In each layer in 
the crystal, the HOMO ((f)) of each molecule splits into 
bonding and antibonding ones due to the overlap with 
its four adjacent inequivalent molecules to form VB1 and 
VB2 Bloch states, respectively. We can then symbolically 
express them, at T, as Vvbi = </>homo - '/'homo anci 
^VB2 = 0homo + ^homo> as confirmed in Fig. 1 (a). 

A B 

We also note that -0VB3 and ^vb4 are mainly <Phomo-i 
in character, so that all the four valence bands mainly 
derive from HOMO and HOMO— 1. If we now turn to 
the conduction bands, they are the 103rd (CB1) - 106th 
(CB4) bands with a width of 0.39 eV, separated from 
the upper ones by 0.15 eV (Fig.l (a)). A similar analysis 
shows that the Bloch states V'cbi - "0CB4 mainly derive 
from 0l{imo(+i) as summa rized in Fig. 2 (a). 

Downfolding — To build a tight-binding model for the 
electronic structure of the pristine picene, we have con- 
structed maximally localized Wannier orbitals (WO's) 15 
for the four conduction bands. The DFT calculations 
were performed with the Tokyo Ab initio Program 
Package (TAPP) 16 and the Troullier-Martins norm- 
conserving pseudopotentials 17 within LDA. The Bloch 
wave functions were expanded in terms of plane waves 
with an energy cutoff of 49 Ry. The features of the cal- 
culated band structure of the conduction bands (Fig. 3 




Fig. 1. (a) Calculated band structure of undopcd crystalline 
picene. The horizontal axis is labeled by (a*/2, b* /2, c*/2), while 
the origin of the energy is set to the conduction band bottom. 
Insets show the highest occupied molecular orbital and the VB1 
and VB2 Bloch wave functions at T. Bottom panel is a blowup of 
the conduction bands (shaded region in the top panel) , with the 
density of states per molecule shown in the solid curve while the 
dashed curves the contributions from the four bands. Horizontal 
lines indicate Ep for x = 1, 2 and 3. (b) The Fermi surfaces when 
E-p is shifted to the situation in K x picene with x = 1,2 and 3 
with rigid bands, (c) The electron density viewed along the long 
molecular axis. 



(b) ), with the atomic positions optimized above, are iden- 
tical to those obtained by QMAS within the numerical 
accuracy. Setting the energy window to cover the four 
conduction bands, we have obtained four maximally lo- 
calized WO's, two of which (wxi,wah) are centered at 
molecule A while the other two (wbu w^h) at molecule B. 
The wave functions w^^WBi (Fig. 3 (a)) have the same 
shape around each molecule for each of i = l,h with 
1(h) denoting the lower (higher) energy states. The band 
structure in Fig. 3(b) in the tight-binding model construc- 
ted from the WO's reproduces accurately the DFT- 
LDA band dispersion. The transfer integrals, Uj(d) = 
(wi(r)\Hjis\wj(r — d)), are defined as the matrix ele- 
ments of the Kohn-Sham Hamiltonian H^s between the 
WO's separated by d — dia + d,2b + d^c. The transfer in- 
tegrals with significant magnitudes are depicted in Fig. 3 

(c) . If we examine the shapes of u>i and Wh, they can be 
represented approximately as </>lumo ± </>lumo+1) with 
Wh having a slightly larger Wannier spread than that of 
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Fig. 2. (a) The lowest unoccupied molecular orbital (LUMO) and 
LUMO+1 for an isolated picene molecule as compared with CB1- 
CB4 Bloch wave functions at T in the crystalline picene. Solid 
lines indicate the main orbital components, (b) 111th and 112th 
Bloch wave functions at T that have amplitudes in the interlayer 
regions. 

to/. As seen in Fig.2 (a), the LUMO (LUMO+1) are an- 
tisymmetric (symmetric) with respect to the plane that 
bisects the long molecular axis, and wi and Wh are nearly 
mirror images of each other, with a slight polarization 
from the center of the molecule. The matrix element of 
H KS between the LUMO and LUMO+1 in the crystal 
is approximately given by |(0lumo|#ks|0lumo+i)| - 
{Eh—ei) /2 = 11.3 meV, where £i = tu(0). This is not neg- 
ligible against the gap of 67.8 meV between the LUMO 
and LUMO+1 calculated for the molecule, which leads to 
the hybridization of the LUMO- and (LUMO+1 )-derived 
bands. This contrasts with pentacene, for which the cal- 
culated gap in the molecule between the LUMO and 
LUMO+1 is as high as 1.28 eV, and the bands derived 
from them do not hybridize. 




Fig. 3. (a) Maximally localized Wannier orbitals with the lower- 
energy W[ and the higher-energy , constructed for the conduc- 
tion bands, (b) Tight-binding fit (curves) as compared with the 
DFT-LDA band structure (dots) for the four conduction bands, 
(c) Arrows indicate the transfer integrals between the in-plane 
molecules whose magnitudes are larger than 25 meV, with el- 
lipses representing the associated WO's, and overlaps represent- 
ing the difference in the c-axis height. 

Doped solid picene — Let us move on to the doped 
solid picene. In contrast to the case of BEDT-TTF or 
TMTSF, carriers are electrons in K^picene. If we assume 
rigid bands and shift the Fermi energy to a position cor- 
responding to x = 3, CB3 and CB4 intersect the Fermi 
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level, and the Fermi surface, displayed in Fig.l (b), shows 
a double surface structure. 18 However, the rigid band as- 
sumption has to be tested. The measured lattice parame- 
ters 8 of the K 2 . 9 picene are a = 8.707, b = 5.912, c = 12.97 
in A with (3 = 92.77°. One conspicuous feature is that 
the c axis shrinks while the a axis expands upon dop- 
ing. This contrasts with the crystalline pentacene, where 
significant expansions of the c axis are known to occur 
when doped with alkali elements 19 as a result of the in- 
tercalants accommodated in the interlayer region. So this 
suggests that, in the case of picene, the dopants are not 
inserted into the interlayer region. 

Another issure of interest is "interlayer states" in GIC. 
There is a body of studies, both experimental 20 and the- 
oretical, 21 that indicate the interlayer states of graphite 
that reside in the space between the carbon layers are 
relevant in the electronic structure in general and super- 
conductivity in particular in GIC, where the interlayer 
states of the pristine graphite, located above the Fermi 
level, strongly hybridize with the valence states of the 
intercalant atoms and the electronic bands responsible 
for the hybridized states descend to intersect the Fermi 
level. In solid picene, we have found that, while there 
exist interlayer states (where we define a layer as the ar- 
ray of molecules in the ah plane) with amplitudes well 
away from the layers, they sit in energy 1.5 eV above 
the conduction band bottom (Fig.2 (b)). If the dopants 
are accommodated within the layer of the solid picene as 
suggested above, the dopant atoms are then expected to 
have little hybridization with the interlayer states of the 
solid picene. 

To determine whether intralayer insertion of K atoms 
is energetically more favorable than interlayer insertion, 
we have prepared initial geometries by locating a K atom 
on each of 3 x 3 grid points on the ab plane of the pristine 
picene, from which we started optimization. The energies 
of the intralayer insertions are found to be about 0.2 eV 
lower in energy per molecule than those of the interlayer 
insertions. We have therefore performed an optimization 
for Kspicene starting from an initial geometry where six 
K atoms are put in between the molecules (four above 
the end benzene rings and two above the center rings) 
within the layer with the lattice parameters set to those 
observed for K^.gpicene. The angle between the molecu- 
lar planes obtained is now obtuse, 114°, which is due to 
the expanded a and shrunk b axes. Analysis of the resul- 
tant structure shows that the smaller volume of the unit 
cell and the presence of the dopants cause a larger dis- 
tortion of the planar shape of the molecules than in the 
pristine case. The electronic band structure is shown in 
Figure 4 (b) (solid lines). They significantly differ from 
those in the pristine picene, so that the rigid band as- 
sumption is invalid. In Fig. 4 (b) we have also displayed 
(dashed lines) the result for a geometry with the K atoms 
removed but the other atoms and the doping level fixed. 
While the two results differ to some extent, which im- 
plies that potassium-states are hybridized in the conduc- 
tion bands, a stronger origin of the non-rigidity against 
doping comes from the difference in the molecular ori- 
entations, which dominates the overlaps of the tt MO's. 
The Bloch states near the Fermi level, however, turn out 



Letter 



4 J. Phys. Soc. Jpn. 



Letter 



Author Name 



to have still the character of the LUMO and LUMO+1, 
while interlayer states were not found near the Fermi 
level. 

The Fermi surface for the doped system in Fig. 4 (c) 
consists of multiple surfaces, which contain warped, in- 
terwoven one-dimcsional (planar) surfaces and a three- 
dimensional (pocket) surface, some of which involve 
transfers between MO's via K atoms. 

To summarize we have obtained the band structure 
of solid picene with curious features. The geometry op- 
timized here is not a decisive one, since there can be 
many metastable geometries, which was also the case 
with solid pentacene, 22 and further investigations should 
be needed. Future problems include how the features of 
the Fermi surface will have implications on mechanisms 
of superconductivity. 
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Fig. 4. (a) Optimized Kspicene geometry with large balls repre- 
senting K atoms, (b) Band structure (solid curves with E = 
corresponding to Ep), along with that for the K-removed geom- 
etry for comparison (dashed curves with E = put equal to the 
position of Ep when we shift it to the situation for x = 3). (c) 
The Fermi surface for the K3 picene. 
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